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ABSTRACT 



We focus on the influence of the Coriolis acceleration on the stochastic excitation of oscillation modes in convective regions of rotating 
stars. Our aim is to estimate the asymmetry between excitation rates of prograde and retrograde modes. We extend the formalism 
derived for obtaining stellar p- and g-mode amplitudes (Samadi & Goupil 2001, Belkacem et al. 2008) to include the effect of the 
Coriolis acceleration. We then study the special case of uniform rotation for slowly rotating stars by performing a perturbative analysis. 
This allows us to consider the cases of the Sun and the CoRoT target HD 49933. We find that, in the subsonic regime, the influence 
of rotation as a direct contribution to mode driving is negligible compared to the Reynolds stress contribution. In slow rotators, the 
indirect effect of the modification of the eigenfunctions on mode excitation is investigated by performing a perturbative analysis of 
the excitation rates. The excitation of solar p modes is found to be affected by rotation with excitation-rate asymmetries between 
prograde and retrograde modes of the order of several percent. Solar low-order g modes are also affected by uniform rotation and their 
excitation-rate asymmetries are found to reach 10 %. The CoRoT target HD 49933 is rotating more rapidly than the Sun (£2/£2 Q ~ 8), 
and we show that the resulting excitation-rate asymmetry is about 10 % for the excitation rates of p modes. We demonstrate that p 
and g mode excitation rates are modified by uniform rotation through the Coriolis acceleration. A study of the effect of differential 
rotation will be presented in a forthcoming paper. 
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1. Introduction 

Internal dynamical processes in stars and wave excitation, prop- 
agation, and induced transport can be strongly influenced by ro- 
tation. Those mechanisms modify stellar internal structure and 
evolution with significant co nsequences for example for galac- 
tic evolution (iMaedej 120094) . The impact of rotation on stars 
is now studied by including models of internal tr ansport pro- 



cesses in stellar evolutionary codes (see for instance | Talon et aTl 
1997tlMaeder & MeynetfeOOOtlEspinosa Lara & Rieutordl2007f 



Decressin et al. I I2009t iMaederl 120091 and references therein). 



Asteroseismology is a lso being increasingly developed with re- 
sults f rom the CoRoT dMichel et al . 2008a b; lAppourchaux et all 



20081) 



120081) and KEPLER (IChristensen-Dalsgaard e~ 
sions, which place constraints on stellar modeling. Those spa- 
tial missions allow us to study stars that are slow as well as very 
rapid rotators. 

Since the pi o neerin g works of lUlrichl (119701) and 
iLeibacher & Steinl (Il97ll) . which led to the identification 
of the solar five-minute oscillations as global acoustic standing 
waves (p modes), the Sun internal structure has been determined 
from the knowledge of its oscillation frequencies. One of the 
remaining key issues is the detection and identificat i on of 

^ 20021: 



gravity modes (|Appourchaux et al 
Turck-Chieze et al.l I2004T 



2000t iGabriel etaL|_ 
Garcia eTalT 120071: iMafhur et al.l 



2007t IGarcfa et al.l l2008al) for determining the rotation profil e 
in the nuclear region (IMafhur et al.ll2Tj0lft IGarcfa et al.ll2008bl) 
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Oscillation modes are indeed crucial for probing the interior of 
rotating stars. 

Stochastic excitation of radial modes by turbulent convec- 
tion has been inves t igated by means of several approaches 
dGoldreich & Keelevl 119771: [Goldreich et alj Il994t iBalmforthl 
ll992tlSamadi & GoupiJ200lHChaplin et al.l2005l) . These meth- 
ods differ from each other in the nature of the assumed exci- 
tation sources, the adopted simplifications and approximations, 
and by the way that the turbulent convectio n is described (see 
reviews bv lStein et al]|2004t iHoud ek 2006). Two major mech- 
anisms have nevertheless been identified as driving the reso- 
nant p modes of the stellar cavity: the first is related to the 
Reynolds stress tensor and, as such, represents a mechanical 
source of excitation; the second is caused by the advection of tur- 
bulent fluctuations of entropy by turbulent motions, and as such 
repre s ents a thermal source of excitation dGoldreich & Keelevl 
mrh ISamadi&Goupiill200l . Samadi & Goupil (2001) pro- 
posed a generalized formalism, taking the Reynolds and entropy 
fluctuation source terms into account. The satisfying agree- 
ment between modelin g and observational data for the Sun 
dBelkacem et al .1 [2006allbh permitted us to go a step further and 
inves tigate the excitation of non-r adial modes in the non-rotating 
case (Belk acem et al J2008L 120091) and now the effect of rotation. 

Our motivation is to investigate the effect of rotation on the 
mode excitation rates rather than the frequencies. We then focus 
on the excitation rates of stochastically excited modes for which 
several issues can be addressed. For example, is the excitation 
rate of a non-axisymmetric mode (m + 0) the same as for an 
axisymmetric one (m = 0)? Are prograde and retrograde modes 
excited in the same manner and what are the consequences? We 
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pay attention to the Coriolis acceleration effects in stars, neglect- 
ing the centrifugal acceleration-induced effects such as star de- 
formation. Our first objective is to determine whether or not uni- 
form rotation can drive the mode efficiently, and our second is 
to evaluate the excitation-rate asymmetry between prograde and 
retrograde modes induced by the perturbation of the eigenfunc- 
tions by uniform rotation. The effect of differential rotation on 
the mode excitation rates will be addressed in a forthcoming pa- 
per. 

The paper is organized as follows. Section 2 introduces the 
general formalism, and a detailed derivation of the Reynolds, en- 
tropy, and rotation-induced source terms is provided. In Sect. 3, 
the formalism is applied to solar spheroidal modes. The special 
case of slow rotators, the Sun, and the CoRoT target HD 49933 
are then investigated, and the results are discussed. Some con- 
clusions are presented in Sect. 4. 



2. Turbulent stochastic excitation 

2.1. The inhomogeneous wave equation 

We derive the inhomogeneous wave equation by taking into 
account the Coriolis acceleration and differential rotation. The 
fluid velocity field (v) is divided into the terms 



v = u + r sin 9 Q(r, 9) e$ . 



(1) 



where Q(r, 9) is the rotational angular frequency assuming an 
axisymmetric rotation, r sin 9 Q.(r, 9) e$ is the velocity field as- 
sociated with rotation, u is the velocity field associated with 
the turbulent convective motion and waves, and (r, 9, (f>) are the 
usual spherical coordinate with their associated unit vector ba- 
sis {£k}k={r,9,6y The rotation axis is chosen so as to coincide with 
the 9-0 axis of the spherical coordinates system of an iner- 
tial frame of reference. In this work, meridional circulation is 
ignored. 

The equation of mass conservation and motion in the pres- 
ence of axisymm etric rotation, can be written as follows (e.g., 
lUnno et al . 1981): 



(2) 



d (pa) 



at 

= pg-VP, 



+ V : (puu) + p 



du 

£1— +2ilxu + rsm6u ■ Vfle,* 

0(f> 



(3) 



where u is the velocity, p is the density, £2 = Q(r, 9) e z is the 
rotation velocity, e z is the unit vector along the rotation axis, g 
is the gravitational field, and P is the pressure. We note that the 
centrifugal force is neglected. 

To go further, all physical quantities are divided into an equi- 
librium one and a perturbation. The subscripts 1 and denote 
Eulerian perturbations and equilibrium quantities, respectively, 
except for velocity where the subscript 1 has been dropped for 
ease of notation. In the following, the velocity field u is divided 
into two contributions, namely the oscillation velocity (v osc ) and 
the turbulent velocity field (it,), such that it = v osc +it,. Then, tak- 
ing the temporal derivative of the equation of motion and using 
the mass conservation equation, one then obtains 



(d 2 \ 

I ^2 ~~ ^ n I V ° sc + ^ osc 



s,, 



(4) 



where Lq is the linear operator that in presence of rotation be- 
comes 

Ln = V [a s v osc ■ Vs + c 2 V (p v OS c)] - ■ (pov osc ) 

„<5 2 v osc <9v osc . dv osc 

~Po Q a , a , - 2 PoQ x ~5 porsmO— — ■ vCle^. (5) 

dtd<p ot at 

The operator C osc involves both turbulent and pulsational ve- 
locities and contributes to the linear dynamical damping (see 
ISamadi & GoupifeOOlL for details): 

d r<9(p,v osc ) „dv osc 
= h: + 2 V : (p()V OS c«r) + p^^— 

+ 2p,n x v osc + p,r sin 9 (v osc • VO) e 

+ V (ff s v osc ■ Vii + c;V • (p,v osc )) ] . (6) 

Finally, the S t operator that contains the source terms of the in- 
homogeneous wave equation (Eq. ©) is given by 



S, = -—V : (p u,u,) + V{a s ii t ■ Vs f ) + Sn + S M 
ot 



and 



p, IQ-^-h, - 2Q x u t — r sin 9u t • V£lc^ 



(p,g\) + V [c]V ■ (p,H,)] - #V ■ (p,u t ) 



at 
a 2 



(7) 



(8) 



(9) 



where g\ is the perturbation of the gravitational acceleration. 
The first two terms of Eq. (0 correspond to the Reynolds stress 
and entropy contributions, respectively. The three following 
terms are contributions associated with rotation. Eventually, as 
shown by Sa madi & Goupill(l2001l) . the terms in Sm do not con- 
tribute significantly to the excitation and are thus neglected as 
well as the linear termfl (£<). 

2.2. Mean square amplitude for uniform rotation 

Using Eq. @, the next step is to determine the mean square 
amplitude of v osc . The procedure is the same as described in 
iBelkacem et al.l (120081) . The wave velo city field is relate d to the 
displacement by means of the relation dUnno et alJI 19891) 



v osc = A [icrij - (£ • V£l)r sm9e^\ e il 



(10) 



where cr — ojq + mCl, and wo is the mode frequency without 
rotation. For uniform rotation, it reduces to 



v osc =Aicrfr ttrt , 



(11) 



where cr is the eigenfrequency, £(r) is the displacement eigen- 
function in absence of turbulence, and A(t) is the amplitude due 
to the turbulent forcing. In the presence of rotation, the wave 
displacement (£) is expressed as 



f m * 



1 dY ( ,„ 
sin 9 d(p 



i A,m 1 0Y Cm c Jn dY^ m \ -i 



(12) 



1 Linear terms are defined as the product of an equilibrium quantity 
and a fluctuating one. 
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where £ r , and £7- are the radial, horizontal, and toroidal com- 
ponents of the displacement eigenfunction, respectively. Note 
that in the following we do not use the upper-scripts f,m on the 
eigenfunction components for ease of notation. Each mode is 
also labelled with a radial order n, which we also omit. 

The power ( P ) injected into each mode with given («, £, m) 
is then related to the mean-squared amplitude ( < \A\ 2 > ) by 



P = n < \A\ 2 > Icr 2 



(13) 



where the operator <> denotes a statistical average performed on 
an infinite number of independent realizations, 77 is the damping 
rate, and / is th e mode inertia- 
Fo llowing ISamadi & Goupill d200ll) and iBelkacem et alj 
one then obtains the mean square amplitude for each 

mode as 



< \A\ 2 >-- 



1 



817(0-/)" 



C c ) 



(14) 



where C\ is the Reynolds stress contribution, C 2 a contains the 
contributions related to the Coriolis acceleration, the Doppler 
term, and one related to the differential rotation, Cj corresponds 
to entropy fluctuation contributions, while C c represents the 
cross-source terms, i.e., the interferences between the different 
source terms. 



2.3. Reynolds stress contribution 

Following the formalism of IBelkacem et alj d2008l) . we develop 
the turbulent Reynolds contribution (see Appendix lAl for a de- 
tailed derivation), which becomes for a given (I, m) 



C 2 R = I6n 4 J drr 2 Pl) R(r)S R (cr), 



(15) 



and 

R(r) = ^ 
v ; 15 



dr 



44 
+ — 

15 



5 I r dr 



■ cc. 



+ 1/ 



ii K+|sP) _|(» + , 

2/ 16 4 8 



- —L 
5 



in 

r 

€t 2 ( 



'—L 2 {l}-2)--T (M --L 2 



dr r 



+ cc 



where 



dfe 1 

+ -(£-&) and 
dr r 



dr r 



(16) 



(17) 



and 



\m\(2(+ 1) 



f,|m| 



[^+l)-(m 2 + l)] (18) 



and L — 1(1 + 1). Furthermore, 



SrQt) = e2 (V J dcj Xk (cj + cr) Xk(o)) , (19) 

where (k, a>) are the wave number and frequency of the turbu- 
lent eddies, and E(k, <d) is the turbulent kinetic energy spectrum, 



which is expressed as the product E{k)xk{<jS) for isotropic turbu- 
lence (ISteinlll967h . A detailed disc ussion of the t empor al corre- 
lation function ixk) is addressed in lSamadi et al.l (12003b . 

Note that in absence of rotation (i.e., D. - 0), the toroidal 
component of the eigenfunction £7- vanishes in Eq. (fT2l and for 
C 2 R and R(r) we recover the ex pressions given by Eqs. (22) and 
(23) of IBelkacem et al.l (120081) . From Eq. £[6j, additional terms 
are found to appear through the toroidal component of the eigen- 
function. All are found to be positive regardless of t and m, im- 
plying an increase in the excitation rates. 

We emphasize that rotation is understood to create 
anisotropics in the Reynolds stress tensor, then off-diagonal 
terms (e.g., iKumar et al.ll 1 995t iMi esch 2005). An adapted spec- 
tral description of turbulent convection including the effect of 
rotation is thus required to compute Eq. ( fT~9b and is beyong the 
scope or our study. 



2.4. Entropy fluctuation contribution 

As sho wn by ISamadi & Goupill (1200 ll) and IBelkacem et al.l 
d2006bl) . the Reynolds stress contribution is not the unique 
source of excitation but one has to account for the excitation 
by the entropy contribution to reproduce the excitation rates for 

solar radial p modes. 

Following IBelkacem et alj (|2008s), the entropy source term 
depends on the mode compressibility that can be estimated as 



I 



r- dr v ' r 



(20) 



where Q is the solid angle, and the spherical harmonics are nor- 
malized following Eq. dA.2) . 

Hence, from Eqs. [T2l and |20l the divergence of the toroidal 
component, which is the curl of the spherical harmonic, van- 
ishes. Consequently, one obtains the same result as for taking 
only the poloidal contribution into account. The final expression 
fo r the contribution of en tropy fluctuations remains the same as 
in lBelkacem et al.l (120081) . i.e., 



C 2 - 



j ' d 3 x a 2 



(Ai + B e )S s ((r), 



(21) 



where "7Y is the anisotropy factor introduced in lSamadi & Goupill 
(|2001|) . which, for the current assumption (isotropic turbulence), 
is equal to 4/3. In addition, 





1 


A t = 


7- 




1 


B ( = 






r 




1 


Dc = 


7- 


where 





D ( 



dr 



d (In I a s |) dD t 



din r 



dlnr 



Ss (0-) = j ^E(k)E s (k) j du X , 



k (cr + u>)xk(w) ■ 



(22) 
(23) 
(24) 

(25) 



In contrast to the Reynolds contribution expression Eq. dToT l, the 
entropy one is not directly modified by rotation. Nevertheless, 
this contribution can be influenced indirectly by means of the 
modification of the radial and horizontal components of the 
eigenfunctions (^,. and %h) by the Coriolis acceleration. 
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2.5. Rotational contributions 



3. Application to spheroidal modes of slow rotators 



The rotational contributions in the inhomogeneous equation 
(Eqs.|4]and[8ll are 

- the contribution related to the Coriolis acceleration 



d d 
- — (2 p,n x u t ) = -2£l x - (p t u t ) , (26) 
ot Ot 



where we neglect the time variations in the angular velocity 
on a dynamical time scale. 
- the contribution related to the Doppler shift 



d ( „du,\ d I du, 



the contribution related to the differential rotation 



(27) 



-^r(p t r sin 6>«, • VQ)e$ = r sin 6» • VQe,*. (28) 
ot ot 

In this paper, we consider only uniform rotation, hence the last 
contribution (Eq. d28| i) vanishes. Nevertheless, all contributions, 
i.e., from Eq. d26l i to Eq. (1281 1. are proportion al to the perturbed 
mass flux p,u,. A dimensional analysis (see ISamadi & Goupill 
1200 ll for details) shows that all those terms then scale as the 
Mach number to the third (Al 3 ). Compared to the Reynolds con- 
tribution, which scales as M 2 , all rotational contributions are 
negligible in the subsonic regime. For the Sun, this conclusion 
remains valid even for the uppermost layers where M ~ 0.3. In 
addition, the rotational velocity appears from Eq. d26l > to Eq. (l28l 
introducing the ratio Q/cr, which is very small for slow rotators. 

Eventually, one obtains (see Eq. ( fT3T >, Eq. (12Tb , and Eq. 
to Eq. <[27j>) 

C\ = 0(M 4 ) » C\ = 0(M b ) » C RS1 = 0(M 5 ) 



»c 2 n = o(m 6 ) o[-\ , 



(29) 



where Crq is the coupled source term associated with the 
Reynolds stress and rotational contributions. Consequently, in 
the following only the Reynolds stress contribution will be con- 
sidered. 



2.6. Final balance 

We have shown in Sects . |2~31 to 12 . 5 1 that in the presence of uni- 
form rotation, the Reynolds term contribution (Eq. (l5[ ) remains 
the most dominant in the subsonic regime. It can be influenced 
by uniform rotation in three ways; 

- The turbulent velocity field can be modified by the Coriolis 
acceleration, hence affecting the Reynolds contribution in 
Eq. (JT3J by the source term (Eq. (|T9b). 

- The toroidal component of the eigenfunction introduces ad- 
ditional terms in Eq. ( [16) . 

- Eventually, the poloidal components of the eigenfunctions 
are modified by the Coriolis acceleration and will influence 
the Reynolds contribution in terms of Eq. (115) . 



As mentioned in Sect. 12.61 the velocity field can be modified by 
the Coriolis acceleration. However, for slow rotators the rotation 
rate does not significantly affect the turbulent field in the upper 
convective region where modes are excited, provided that the 
ratio of the convective frequency to the rotation rate is higher 
than unity. For the Sun, this requirement is fulfilled in the en- 
tire convective region except in the deepest layers, near the in- 
terface with the radiative region. Nevertheless, the contribution 
of these deep layers do not contribute significantly to the ex- 
citation rates for the modes considered here, i.e., low-order g 
modes and p modes. Hence, in the following we assume that 
the turbulent field, and its spectral dependence, are not affected 
by uniform rotation. Note, however, that for lower frequencies, 
and especially asymptotic gravity-modes, this approximation is 
no longer valid since a significant contribution t o the mode exci- 
tation comes from the deeper convective layers dBelkacem et al.l 
120091) . 

We then consider the effect of the perturbation of the mode 
excitation rates by the Coriolis acceleration associated with the 
modification of the eigenfunctions. In this framework, we use a 
perturbative approach, which is valid for slow rotators and par- 
ticularly for the Sun since we restrict our investigation to rather 
high-frequency p and g modes. The ratio of the mode frequency 
to the rotation rate is still higher than unity. In the Sun, for a typ- 
ical p mode at v = 3mHz, one has u>o/2Q. « 3 x 10 3 » 1 in 
the convective region where modes are excited, and for a solar g 
mode at v — 100 /iHz, this ratio remains high at o>o /2CI « 100 » 
1 . This allows us to use a perturbative approach. 

3.1. Perturbation of the mode excitation rates 

Our aim is to derive an analytical expression of excitation rates 
influenced by uniform rotation using a perturbative me thod. 
Following the classical method given in lUnno et al.l d 19891) . we 
develop the eigenfunction for a given n, t, m to first order 



n +n 



(30) 
(31) 



and for the toroidal part 

- {t + 1) D MM [£l M + ({ + 2)4°^ +1 ] } , (32) 



where 



D 



t,m 



c 2 



i 



t 2 V 4-e 2 - l 



(33) 



The expression for C„' , hC is given in AppendixE and g® v 

are the radial and horizontal components of the eigenfunction in 

absence of rotation. 

As pointed out by iDziembowski & Goodd dl992l) . the con- 
vergence properties of the sum involved in Eq. (IB. 1 b and 
Eq. ( IB. 2b are unclear. As shown by Eq. (IB ,4b . this is particu- 
larly problematic for a dense spectrum such as high-order grav- 
ity modes in the S un. We nevertheless use it for c onvenience. An 
alternative exists (Dziembowski & Goode 1992), which consists 
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0.001 



0.1 1.0 
V [mHz] 

Fig. 1. Top: Excitation rates for the mode I = 1 and m = 
{-1,0, 1), computed using Eqs. ( TT3"|l. ([ T4 |, and (|T5t and using 
the same solar model as detailed in Belkace m et al] (12008b . The 
depression in the excitation rates at v ss 300 /^Hz is associated 
with the mixed nature of the modes that consequently produce 
a higher inertia thus a lower P. Bottom: Bias between prograde 
and retrograde modes, defined in Eq. d37b . for the same modes 
as for the figure in the top panel. 



of computing a modified eigenvalue problem. However, this sec- 
ond possibility makes it more difficult to identify contributions 
to the excitation rates. 

Inserting the decomposition (Eqs. IB.ll to IB. 31 ) into Eqs. [TBI 
and[TJj] we obtain (see AppendixlBlfor the detailed calculation) 

'-'- + -(i?)'S- m 

where 

P(0) = ^ / dmR^S R (u ), (35) 

P M = W J dmBf&SAo*), 06) 

and and corresponds to the perturbative expansion of R 
(Eq. ( fT6] l) given in AppendixlBl Note that the zeroth-order terms 
(°) correspond to the case without rotation. Only the first order 
in (2Q/<x) is considered. Accordingly, the contributions of in 
Eq. (TT3T > are neglected because they are of second order. 



Fig. 2. Bias between prograde and retrograde modes, defined in 
Eq. d37l i. for four values of the azimuthal order m and the radial 
order n — 5 as a function of the angular degree (. The computa- 
tions are performed in the same manner as for Fig.Q] 



We now define the excitation rates asymmetry, between pro- 
grade and retrograde modes to first order such as 



(-) 


( p(l)\ 
N 


Wo/ 


p(0) 



(37) 



From Eq. d37l i. two factors contribute to the asymmetry namely 
the ratio (2Q/w ) and Pf,J,[/P (0) . They are discussed in the fol- 
lowing sections. 

3.2. Application to slow rotators 
3.2.1. The Sun 

Using the same nu merical computation as described in 
iBelkacem et all (I2008I) . we apply this formalism (Eq. (O - [36b 
to the solar case. FigureQ]displays the mode excitation rates for 
the t = 1, p and g modes. It also presents the ratio SP m /P m , 
defined in Eq. fl37| l, which emphasizes the effect of the mode 
excitation rate asymmetry between the prograde and retrograde 
modes. 

It turns out that the excitation rates of acoustic modes are 
modified with an excitation rate asymmetry of the order of the 
percent, which increases toward g modes. We find that the vari- 
ation in the mode excitation-rate asymmetry with frequency is 
caused by the term 2£2/a>o in Eq. (l3"7T l. while the ratio P\2/P^ 

remains of the order of the value of one. Pj^j is dominated by the 
first term in Eq. (|B.12| i, which corresponds to the contribution of 
the radial component of the eigenfunction, for p modes. For g 
modes, the horizontal component of the eigenfunction is also of 
importance and contributes significantly to p| m |. 

For higher values of the angular degree ({), as shown by 
Fig. |2] there are two effects. First, the higher the azimuthal or- 
der m, the higher the mode excitation rates asymmetry, at fixed 
i. This is explained by the perturbation of the mode excitation 
rates being proportional to m in Eq. d34l l. Second, at fixed m, 
the higher the angular degree, the lower the mode excitation- 
rate asymmetry. This behavior comes from the frequency shift 
of high-f modes, since at fixed radial order, the higher the angu- 
lar degree the higher the mode frequency. Hence, for the same 



6 



K. Belkacem et al.: Mode excitation by turbulent convection in rotating stars. I. 




0.001 



0.5 1.0 1.5 2.0 2.5 

V [mHz] 

Fig. 3. Bias between prograde and retrograde modes, for the i = 
1 modes of the star HD49933 (P n =3.4 d) and the Sun (P n =2S 
d). 



radial order the ratio 2Q/wo will decrease with the angular de- 
gree explaining the behavior in Fig. [2] 

3.2.2. The CoRoT target HD 49933 

We now consider more rapid rotators, such as HD 49933. This 
is an F5 V main-sequence star observed twice by the CoRoT 
missiorQ, first during 62 days and more recently for more than 
150 days. The unpre cedented photometric precision achieved by 
the CoRoT mission dMichel et al.ll2008al: lAuvergne et al.ll2009l) 
makes this star a good candidate for the detection of mode 
excitation-rate asymmetry, which requires, as previously men- 
tioned, accurate measurements. This star, indeed, exhibits a sur- 
face rotation period that is shorter th an that of the Sun, Pq « 3. 4 
days (i.e., Q/C2 w 8) as shown by lAppourchaux et al. (2008). 
but still slow enough to ensure that the perturbative approach is 
valid. 

In Fig. [3] we present the same ratio as in Fig.|2]for the I = 1, 
p modes using a mod el of HD 49933 that match e s the seismic 
const raints derived by lAppourchaux et al] d2008l) dGoupil et alj 
120091) . The asymmetry between the excitation rates of m — 1 and 
m = -1 modes is found to reach up to 10%. In terms of mode 
excitation-rate asymmetry, the differences between the Sun and 
HD 49933 is due to a higher ratio (2£l/<oo) in Eq. d37l 

This demonstrates that an asymmetry in terms of mode ex- 
citation rates is more likely to be observable for more rapid ro- 
tators than the Sun, even if, in contrast to the Sun, only low-^ 
modes are observed. 

4. Conclusion and perspectives 

4.1. Conclusion 

We have derived a formalism that models the stochastic exci- 
tation of oscillation modes by convective motions in uniformly 
rotating stellar regions. We have shown that the driving terms, 
due to rotation, that appear in the inhomogeneous wave equation 
are negligible with respect to the Reynolds stress contribution. 

2 The CoRoT space mission, launched on December 27th 2006, has 
been developed and is operated by CNES, with the contribution of 
Austria, Belgium, Brazil , ES A (RSSD and Science Program), Germany 
and Spain. 



We demonstrate that the dominant contribution to the excitation 
rates then comes from the modification of the eigenfunction by 
the Coriolis acceleration. 

The formalism is then applied to low-order g modes and p 
modes of slow rotators, and in particular the Sun and the CoRoT 
target HD 49933. For the Sun, a bias between pro- and retrograde 
waves is found in the excitation rates. For g modes, this bias 
can reach values of up to 10 %. For low-{ p modes, this bias 
is found to be of the order of a percent. The detection of the 
mode excitation-rate asymmetry of individual p modes is not yet 
possible since the observatio nal error bars obtained, f or instance 
by GOLF, are around 20 % dBelkacem et al.ll2006bl) . while we 
search for a physical effect of only several percent. 

For more rapid rotators, such as HD49933, we find that the 
excitation-rate asymmetry of low-^ (I - 1) p modes can reach up 
to 10 %. However, this value is achieved at low frequency, where 
seismic measurements are generally dominated by the granula- 
tion background. In the case of the CoRoT target HD49933, de- 
tection of acoustic modes is limited to the frequency domain v e 
[1.2; 2.5] mHz. In this frequency domain, the asymmetry in P is 
no greater than 2%. In contrast, the 1-cr uncertainties associated 
with P are in the range 30% - 80%, depending on the frequency, 
for the observations completed during the C oRoT initial run 
dAppourchaux et alJl200S ISamadi et al.ll2009l) . For the second 
set of observations of HD 49933 by CoRoT, the 1-cr uncertain- 
ties associated with P is expected to be lower, i.e., in the range 
20% - 40 % (Benomar, priva te communication). Furth ermore, 
current seismic analyses (e.g.. lAppourchaux et alj|2008h do not 
reproduce individual mode multiplets but assume a fixed ampli- 
tude ratio of the different mode multiplets or even assume a fixed 
amplitude ratio of the different I degree. Therefore, despite the 
high precision of the CoRoT instrument, it is presently not pos- 
sible to constraint P for an individual mode multiplet (i.e., for a 
given value of 6 and m). Concerning the Kepler instrument, its 
performance in terms of photon noise level is expected to be a 
factor five lower in terms of power compar ed to that achieved fo r 
the brightest stars of the CoRoT mission dChaplin et al. I l2008h . 
On the other hand, Kepler will observe the seismic targets over 
a much longer period (around 4 years) than the CoRoT mission, 
which we hope will permit us to constrain individual mode mul- 
tiplets. Another way to proceed is to consider a sum of P for a 
given m so as to reduce the actual observational errorbars in both 
the Solar case and that of HD49933. 

We note that mode amplitude is a balance between driving 
and damping. Therefore, asymmetries in mode amplitudes can- 
not be inferred only from excitation rates since some possible 
asymmetry in the mode damping rate can arise. This is not in- 
vestigated here but left to future work. 

4.2. Perspectives 

The effect of uniform rotation on the mode amplitude excitation 
rates presented here is exploratory work that requires further in- 
vestigation and theoretical developments. 

Stellar convection zones are differentially rotating. 
Therefore, the next step would be to take the differential 
rotation into account in both the radial and the latitudinal 
directions. In contrast to uniform rotation, a consideration at 
Eq. dTDb permits us to understand that the driving source terms 
in the inhomogeneous wave equation Eq. (0]i are modified by 
differential rotation. The effect of differential rotation on mode 
excitation rates is the scope of an upcoming paper. 

The regime of rapid rotation should also be addressed. The 
formalism must be adapted to the specific geometry of those 
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stars, since spherical coordinates become inappropriate and one 
may have to take the star deformation into account. In ad- 
dition, the eigenfunctions and frequencies have to be derived 
from an adapted non-perturbative method (e.g., iReese et ail 
120061) since rotation can strongly modify both the eigenfunctions 
and the stellar structure of a star (Rieutord & Valdettarolll997t 
iDintrans & Rieutordl l2000t IReese et al.l 120061) . Furthermore, 
in such a regime new types of waves appear that deserve 
a furt her study, such as inertial waves and gravito-inertial 
waves (iRieutord & Valdettard[l99^ iDintrans & Rieutordll2000t 
iMathis et alj|2008l) . In addition, in contrast to slow rotators, the 
turbulent field is also affected by rotation and the spectral de- 
scription of turbulence must be taken it into account. The as- 
sumption of isotropic turbulence is should then be excluded 
and the spectral properties of the turbulent field be specified. 
Numerical simulations can be of some help. As do ne for solar g 
modes using the ASH code dBelkacem et al. 2009), it is possible 
to assess the turbulent properties of these rotators from numeri- 
cal simulations. 
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Appendix A: Detailed expressions for the Reynolds 
source term 



The eigenfunctions (£) are developed in spherical coordinates 
(e r ,eg,e,f,) and expanded as a sum over spherical harmonics. 
Hence the fluid displacement eigenfunction for a mode with 
given {, m is written as 



1 dY 0n 
sin d(f) 



80 
80 



,T;t,m 



1 8Y tj , 
sin d(p 



where the spherical harmonics (Y( m (6, <p)) are normalized ac- 
cording to 



dQ. 
4n~ 



Y e,m Yf - 1 



(A.2) 



with Q being the solid angle (dQ = sin#d#d0). 



The Reynolds str ess contribution can be written as (see 
iBelkacem etai1l2008l for details) 



C 2 = 7T 2 J" d 3 x (p 2 b*jb lm ) j d'k j dco 



(T iilm +r^)^- X k{cT + uS)Xk{oS) 



k 4 



where 
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(A.4) 
(A.5) 



where the double dot denotes the tensor product. 



To compute the coefficients bg in Eq. (|A.3b . we follow the 
procedure derived by IBelkacem et all (120081) . which infers that 
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The contribution of the Reynolds stress can thus be written as: 

3 j dm jdk jdcj R(r, k) 



C 2 = 47T 



E (k) 
k 



(A.6) 



where we have defined dm = Anr podr. Using the Einstein sum- 
mation convention 



(A . 3) R ir , k) = J^J^ bhblm (T^ + T^). 
Because T l ' lm = T-'"' 1 ', it is easy to show that 

R{r > k) =i d ^ rsr B >J B >"< ( r ° h,+r *o. 

where By = (l/2)(b u + b ft ). 

Using the expression Eq. jAAj for T l ' lm , we write 
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where £ <0) , , ; are so lutions 
out rotation dUnnoet al.ll 19891) . 



We assume isotropic turbulence, hence the k components where £ <0) , , gy, ; are so lutions of the oscillation equation with- 
satisfy 



r kk, r k 2 r 
J = ««J dQ *F 

where dg is the Kro necker symbol for i,j = r,9,<f>. As in C„',,,,z 



LO- 



BelkacemetalJ d2008l) . we then obtain 



(B.4) 
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and 
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Usin g Eq. ()A.6|> to comput e Eq. dA.8b . with a = 1 /3 and /? = 
1/5 (see lBelkacem et ai1l200 8, for details), yields 
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The Coriolis corrective terms are of the order of 2Q./a>o, which 
is here assumed to be small. Inserting Eqs. (lB.UB.2HB.3b into 
Eq. (fT5l ), we obtain 
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where we have defined 
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where 
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and /? (0) and correspond to the perturbative expansion of /? 
given in Eq. ( fT6b where 



Appendix B: First-order perturbation of the 
excitation rates 

We recall the main results about the first-order perturbation of 
a spheroidal mode £ (0) due to the Coriolis acceleration (of fre- 
quency ojq) and establish the perturbation of the excitation rates. 
Following the classical method given in Unno et al. (1989), we 
obtain 
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